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We present two different theories for Raman scattering and Resonant Inelastic X-ray Scattering 
(RIXS) in the low temperature ferromagnetic phase of YTiC>3 and compare this to the available 
experimental data. For description of the orbital ground-state and orbital excitations, we consider 
two models corresponding to two theoretical limits: one where the t2 g orbitals are degenerate, and the 
other where strong lattice distortions split them. In the former model the orbitals interact through 
superexchange. The resulting superexchange Hamiltonian yields an orbitally ordered ground state 
with collective orbital excitations on top of it - the orbitons. In the orbital-lattice model, on the 
other hand, distortions lead to local dd-transitions between crystal field levels. Correspondingly, 
the orbital response functions that determine Raman and RIXS lineshapes and intensities are of 
cooperative or single-ion character. We find that the superexchange model yields theoretical Raman 
and RIXS spectra that fit very well to the experimental data. 
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I. INTRODUCTION 

In many 3d transition metal compounds, like the 
titanates and colossal magnetoresistance manganitcs, 
the low-energy orbital degrees of freedom are (approx- 
imately) degenerate. This leads to all kinds of interest- 
ing phenomena like cooperative Jahn- Teller (JT) distor- 
tions, orbital frustration and strong spin-orbit coupling. 
Orbitals on neighboring orbitally active ions can be cou- 
pled via JT-distortions and via superexchange^^]. Both 
couplings can give rise to an orbitally ordered ground 
state. However, the nature of the orbital excitations on 
top of the orbitally ordered ground state is very differ- 
ent depending on the coupling mechanism. In a pure 
superexchange system, collective excitations called 'or- 
bitons' emerge. These collective orbital excitations are 
coherent waves of excited orbitals, similar to the spin 
wave excitations in spin systems. When this system is 
coupled to the lattice, as is the case when JT-distortions 
are present, interaction with phonons destroys the or- 
biton coherence. The distortions lift the orbital degener- 
acy and the excitations become local dd-excitations. 

The titanates, with a pseudo-cubic perovskite lat- 
tice structure, are good candidates to support orbitons. 
The Ti ions with their 3d 1 configuration have one elec- 
tron in one of the three nearly degenerate t2 g orbitals. 
Since these orbitals are directed away from the neigh- 
boring oxygen ions, the coupling to the lattice is ex- 
pected to be small. Further, it has been showrP^O that a 
superexchange-only model explains many of the ground 
state properties of YTiC>3. Also, there is experimental ev- 
idence that LaTiC>3 is a (superexchange-driven) orbital 
liquicpHU Q n ^hc othgj. hand, local crystal field mod- 
els also well reproduce some of the physical properties of 



the titanates^^l. Both models have their shortcomings 
as well: a Jahn- Teller dominated description is not able 
to reproduce the spin wave spectrum, which is nearly 
isotropic in both spin and real space, while the superex- 
change model has difficulties explaining the experimen- 
tally observed orbital polarizatioEpSEH Consequently, it 
still remains controversial which mechanism dominates 
the orbitals in titanates^. 

In order to resolve this controversy, it is of cru- 
cial importance to compare recent Raman and Reso- 
nant Inelastic X-ray Scattering (RIXS) experiments on 
titanates^K^U to both of the competing theories. 

The experimental observation of orbitons is a difficult 
task, as is underscored by the heavily debated Raman 
measurements on LaMnO^H^H. Because the hallmark 
of collective excitations is dispersion, a much better tech- 
nique to directly probe orbitons is RIXS. With its energy 
resolution drastically improved over the last few years, 
RIXS now offers a whole new way of accessing the ele- 
mentary excitations of solids, complementary to for in- 
stance Raman and neutron scattering. In this paper we 
analyze recent Raman and RIXS spectra^^U f or YTiC>3 
from the point of view of a superexchange-only model 
and the alternative extreme of a completely local, lat- 
tice distortion dominated model. We find that while the 
orbital-lattice model can be finetuned to capture some 
aspects of the observed spectra, the collective superex- 
change model yields a much better overall description of 
the Raman and RIXS data. 

The paper is organized as follows: Sec. [IT] compactly 
reviews previous work on YTiC>3 and introduces the su- 
perexchange formalism and the local crystal field model. 
Sections |III| and |IV| deal with the theory of Raman scat- 
tering and RIXS respectively, in both the superexchange 
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and crystal field models. 

II. TWO MODELS OF YTIO3 

For the existence of collective excitations of or- 
bitals, the so-called orbitons, it makes a difference 
whether the orbital order is driven by JT distortions or 
superexchange^!. For large JT distortions, the crystal 
field splitting is large and a local picture applies: the 
collective nature of the orbital excitations characteristic 
of orbitons is lost. In materials where the orbital-lattice 
coupling is small, the superexchange interactions between 
orbitals can dominate over crystal field splittings due to 
lattice distortions. The Ti ions have a 3d 1 configura- 
tion, and the octahedral crystal field induces a splitting 
between the higher energy e g and lower energy t^g lev- 
els. Because the t2 g orbitals are not directed towards 
neighboring oxygen ions, they are not expected to couple 
strongly to lattice distortions. 

Building on this assumption, one can derive a superex- 
change Hamiltonian starting from a Hubbard model. Be- 
low, we follow Refs. [Hand [5] closely. By symmetry, the 
hopping term connects, for instance, the zx to zx and yz 
to yz orbitals along the z-direction (c-axis) via the inter- 
mediate oxygen 2p^ states, xy orbitals are not coupled 
along this direction. In the limit of large on-site Coulomb 
repulsion U, this leads to a superexchange interaction 
that depends on the spatial direction of a bond, and the 
resulting model is intrinsically frustrated: on any given 
ion, there is no orbital that minimizes the bond energy 
in all directions simultaneously. 

Because YTi03 is ferromagnetic at low temperature 2 -^ 
(T c w 30 K), we restrict ourselves to the completely fer- 
romagnetic part of the Hilbert space. Then one obtains 
the simple Hamiltonian 

m v J 

with the orbital exchange integral J or b = ^i^S-Bi where 
r*i = 1/(1 — 3Jh/U) ~ 1.56 parametrizing Hund's rule 
coupling and J$e = At 2 /U is the superexchange constant 
derived from the Hubbard model. The operator 
depends on the direction 7 of the bond ij. For example, 
in the z-direction we have 

A[f = n ati n a ,j + n bii n bJ + a\b l b^a j + b\a i a^ j b j . (2) 

The operators a' , w and create an electron in the yz- 
, zx- and xy-orbital, respectively, and n a = a^a. The 
Hamiltonian can also be written in terms of interacting 
effective angular momenta 1 = 1, operating on the t 2g 
triplet. Because of the orbital frustration, these can form 
a myriad of different classical ground states. Refs. [U |S] 
conclude that a 4-sublattice quadrupole ordered state is 
favored, in which the orbitals 

H'c)= -^(\d yz )±\d zx }±\d xy )) (3) 



are condensed. The signs ± alternate between the sub- 
lattices, such that nearest-neighbor orbitals are orthogo- 
nal, supporting ferromagnetic order. On top of this con- 
densate, two species of orbitons can be created, loosely 
speaking by populating cither one of the two orbitals or- 
thogonal to tj) c . The orbiton spectrum has 3N 1 / 3 Gold- 
stone modes (where N is the total number of Ti ions), 
because the number of orbitals of a specific "color" is 
conserved in the plane in which it is lying. However, 
in YTi03 the TiOg octahedra are tilted. Because of 
this, hopping between different t 2g orbitals is now no 
longer symmetry forbidden, and the conservation of or- 
bital "color" is violated, removing the Goldstone modes. 
When also some anharmonic terms of the Hamiltonian 
are taken into account on a mean field level, the orbiton 
dispersion becomes^ 

Wl/2,k = ^ZjJ mh {\ - (1 - 2e)(l - 2/)( 7l , k ± Kk ) 2 

-2(e-/)( 7ljk ±K k )} 1/2 , (4) 

where we use the signs + and — for wi ik and W2,k re- 
spectively. Further, ^fZ e Z f w 1.96, / « 0.086, e w 0.18, 

71. k = (c x + Cy + c 2 )/3 and K k = y' 7 2 k + 7 2 k with 

72. k = V3(c y - c x )/6 and 7 3:k = (2c z - c x - c y )/6 
with c a = cos k a . Eq. Q describes the collective orbital 
modes that disperse up to energies of 2 J or b and have a 
gap of approximately J or b ■ 

In the second orbital model for YTiC>3 that we con- 
sider, lattice distortions domin ate over superexchange in- 
teractions. Pavarini et aZpE] did a DMFT+LDA calcu- 
lation and found that lattice distortions of the GdFeC>3- 
type lift the orbital degeneracy. They also obtained four 
sublattices. The resulting local eigenstates of the t 2g sys- 
tem areP 

|1) = 0.781 \yz) - 0.073 \zx) + 0.620 \xy) (5) 
|2) = -0.571 \yz) + 0.319 \zx) + 0.757 \xy) (6) 
|3) = 0.253 \yz) + 0.945 \zx) - 0.207 \xy) (7) 

for sublattice 1, with corresponding orbital energies t\ = 
289 meV, e 2 = 488 meV and e 3 = 620 meV. This 
yields excitation energies lo\ = e 2 — (\ = 199 meV, 
^2 = £3 — ei = 331 meV. The orbital states on the 
other sublattices can be obtained from lattice symmetry 
considerations^. Superexchange processes are treated 
as a perturbation in this model, broadening the states 
generated by lattice distortions. This pict ure is also sup- 
ported by other theoretical workP2I2SM] 

It is possible to rotate the axes on each of the sub- 
lattices in such a way that in the new coordinates, the 
eigenstates are still given by Eqs. ([5| through Q: 



subl. 1 : (x, y, z) i-> (x, y, z) (8) 

subl. 2 : (x,y,z) h-> (y,x,z) (9) 

subl. 3 : (x, y, z) i-> (x, y, -z) (10) 

subl. 4 : (x,y, z) i-> (y,x,-z). (11) 
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Correspondingly, the orbiton operators transform as fol- 
lows: 



subl. 1 
subl. 2 
subl. 3 
subl. 4 



(a,b,c) i-> (a,b,c) 
(a,b,c) i-> (b,a,c) 
(a,b,c) i-> (-a,-b,c) 
(a,b,c) h-> (-6,-o,c). 



(12) 
(13) 
(14) 
(15) 



III. RAMAN SCATTERING 

In the search for orbitons, Raman scattering has been 
an important tool for experimentalists. After the con- 
troversial first observation of orbitons in LaMnC>3 24 26 , 
the titanates now seem to be a more promising candi- 
date. In addition to the reasons mentioned in previous 
sections, recent Raman data by Ulrich et alW^ should be 
noted, which shows a striking temperature dependence: 
the spectral weight of the 235 meV peak in YTiC>3 in- 
creases dramatically when temperature is lowered. This 
can be naturally explained by collective orbitons: as tem- 
perature drops, the orbitons gain coherence and the spec- 
tral weight increases, analogous to two-magnon Raman 
scattering in the cuprateiP^. From the local dci-excitation 
point of view, temperature should not affect the inten- 
sity of local transitions between crystal field levels. Also, 
Ulrich et al. found that the polarization dependence of 
the spectra is hard to reconcile with the local excitation 
picture a result that we will reproduce below. In opti- 
cal data^, a peak is seen at the same energy and was 
ascribed to orbital excitations. 

Earlier theoretical work on Raman scattering in the 
titanates^ built on the assumption that JT-distortions 
determine the symmetry of the orbital order. In this 
paper, we investigate the Raman spectrum of YT1O3 
in both the lattice distortion and superexchange frame- 
works laid out in Sec. |TTJ We start out wit h the Loudon- 
Fleury effective Raman scattering operatoi ! 32 * 33 ! 



A 



(7) 



7) 1 



(16) 



where the usual spin exchange Hamiltonian has been re- 
placed by the orbital Hamiltonian of Eq. Q. €ij are the 
electric field vectors of the in- and out-going light, 6ij 
connects nearest neighbors i and j. The physical picture 
is that the light induces an electric dipole transition to 
the intermediate state where a 3d ti g electron ends up 
on a neighboring Ti ion, after which one of the electrons 
of this now doubly occupied site can hop back in another 
transition. In this process, the two involved electrons can 
end up in different orbitals, resulting in a two-orbiton ex- 
citation, in full analogy with two-magnon Raman scat- 
tering in the cuprates. As the light forces the electrons to 
perform a superexchange process independently of the in- 
trinsic coupling mechanism of the orbitals, this effective 



Raman operator holds for the lattice distortion model 
too. 

With this scattering operator, we calculate the Raman 
spectrum for the superexchange model. Similar calcu- 
lations have been done before in the context of Raman 
scattering on orbital excitations in vanadates^. Adopt- 
ing the geometry used in the experiment of Ref. |21j . 
we take the electric field vectors to be in the plane 
parallel to the [110]- and [001]-directions: e i(f) 

sin -^sinflj^cosfy/)) where 6» i(/) is the angle 
the electric field vector makes with the c-axis. Through- 
out this paper we use a coordinate system in which the 
nearest neighbor Ti-Ti bonds are parallel to the coor- 
dinate axes. Substituting into Eq. ( 16 1 and using that 



7i 2 



is a conserved quantity in the superexchange 



model and that H 
man scattering 



oc 



|0), we find for inelastic Ra- 



/»' \ ( cib ti, cos 6f — - sin tii sin tif I 



(17) 



where the sum is over bonds in the c-direction only. Per- 
forming the transformations mentioned in Sec. [TTJ con- 
densing ip c and Fourier transforming, we obtain 

E i(c) - 3 E H - ^ ^ + f (a * - 5 * )x 



(hi). 



(«- k - bl k ) + -f (a_ k - 6_ k )(a k - b k )\ (18) 



where only quadratic terms in the operators are retained. 
Linear terms do not appear. Next, this result is Bogoli- 
ubov transformed according to 

a k = w k ch 0i M a 1M + w k ch 6> 2 ,ka 2 , k 

- w k sh 0i,ka^_ k - w k sh &2,k a l,~k> ( 19 ) 

bit = -w k ch ^i,kQ!i )k + "kch 6> 2: ka 2 , k 

+ i; k sh 6>i ik ix[ ; _ k - u k sh 02,kor 2 _ k , (20) 

where the indices 1,2 refer to the orbiton branch. This 
transformation diagonalizes Hq up to quadratic order if 



u k 



72, k 

2K k 



w k = sign(7 3ik ) 
tanh20i( 2 ), k = 71, k ± K k . 



72, k 

2K k 



(21) 

(22) 
(23) 



The effective Raman scattering operator now cither pro- 
duces two orbitons or scatters single orbitons already 
present in the initial state. At zero temperature, the 
initial state has no orbitons (in "linear orbital wave the- 
ory", i.e. if we neglect orbiton-orbiton interactions), so 
we keep only the two-orbiton creation part of y\- ^ A^ 
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in Eq. ( 18 I 



I [K u + w ) 2 ( c - ch 20 i - sh 2 ^)} «U«l-k 

k 

+ {(« - v) 2 (c z ch 26*2 - sh 26> 2 )} 4,k a 2,-k 

+ 2 {(« 2 - i; 2 )[sh {e l + e 2 ) - Cz ch {e 1 + e 2 )]}a\ M a\ ^ 

(24)" 

where c z = cos & z and the index k is implied on every 
u, v, 9i and 9 2 - 

The cross section at zero temperature now is 



d 2 a 
dcudfl 



f 



(25) 



with / labelling the two-orbiton final states with energy 
u)f. The corresponding matrix elements are given by 
Eq. (|24|). 



Because there are orbiton-orbiton interaction terms in 
the Hamiltonian which are neglected in "linear orbital 
wave theory" , we introduce a phenomenological orbiton 
damping of 7 = 30 meV. Also, broadening from other 
sources such as interaction with phonons and magnons 
can be mimicked this way. 

The result is displayed in Fig. [T] compared to the 
data from Ref. |21j . In the superexchange model, only 
two-orbiton creation processes contribute to the Raman 
spectrum. The best fit is obtained for J or b = 65 meV, 
close to the value estimated in Ref. [5] from magnon 
data of YTiOgEO. Including orbiton-orbiton interactions 
will probably reduce the peak energy (in analogy to two- 
magnon Raman scattering), increasing the fit parameter 
Jse- 

The local model of YTiC>3 also yields Raman spectra 



via Eq. (16 1. In this model, the orbital order makes the 
c-direction different from the a and b ones. Therefore, all 
bond directions are considered separately. For technical 
convenience, the rotations Eqs. ( 12 1 through ( 15 1 are first 



performed. Bonds in the c-direction connect sublattice 1 
to sublattice 3, and 2 to 4. Both these bonds give the 
same contribution to the Raman operator: 



(26) 



^ (na^riaj + ribjTibj + <4Mj a j + h ',","', h . 



+ 2 ( n a,i + n b, 3 ) ) ■ 



Note that the expression is symmetric in i,j. Similarly, 
for the a- and 6-directions, we obtain again the same 
contribution for both bonds with i 6 sublattice 1 and 



j £ sublattice 2, and for bonds with i G 3 and j E 4: 

Y ( n hi n a,j + n c ^n c ,j + b\c i c\a j + ct^atc^ 



(id) a 



+ ( n a,i + n b,j) ) J 



(n a ,inb,j + n Ct in c j + atc^bj + cja^jc^ 



(ij)„ 



+ 2 ( n b,i + n a,j) 



(27) 



(28) 



In general, these operators give rise to final states with 
one and two c?d-excitations. Using the local wave func- 
tions proposed in Ref. [28], final states with one dd- 
excitation cannot be reached in (z, z) polarization config- 
uration, in agreement with the findings of Ref. |29j . Be- 
cause the wave functions Eqs. ^ through Q of Pavarini 
et al. are close to these states, there is little single dd- 
excitation weight (in particular in (z,z) polarization), 
and the spectrum is dominated by double dci-excitations. 
In the numerical calculations of the Raman spectra, the 
same broadening of 7 = 30 meV as above is included. 

The resulting Raman spectra are shown in Fig. [T] to- 
gether with the experimental data. The experimental 
data peaks around 230 meV in the (z, z) polarization 
configuration shown here. In the experiment, other con- 
figurations give very similar line shapes, with the maxi- 
mum shifting around no more than ~ 40 meV. The in- 
tensity is strongest when both in- and outgoing polariza- 
tions are directed along one of the cubic axes 2 ^, i.e., in 
the zz, xx, yy polarization geometries. 

Even though we have included possible orbiton-orbiton 
interactions only as a phenomenological damping, the su- 
perexchange model gives a very good fit to the exper- 
imental line shape: it reproduces a single peak without 
internal structure at approximately the right energy. The 
cubic isotropy of the superexchange model is in agree- 
ment with experiment, as noted in Ref. |21j . 

An interpretation of the Raman spectrum in terms of 
local crystal field excitations is problematic. Not only 
is the predicted strong polarization dependence of the 
intensity (a stark contrast between the c-axis and the 
a, 6-axes) opposite of what is seen in experiment (which 
obeys cubic symmetry^), the suppression of the single 
dii-excitations with respect to double excitations leads 
to a wrong prediction of the peak energy. We tried to 
include corrections to the Raman operator from nondi- 
agonal hoppings between i 2 g orbitals but this did not im- 
prove the fit. Also, to blur the multiple peaks together 
into one peak, a large broadening is needed. Finally, 
the temperature dependence of the peak as observed in 
Ref. [5T] is difficult to explain in the context of local dd- 
excitations. 
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Raman spectra 

SE^ (x,x) and (z,z) 
Local, (z,z) 
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0.2 0.3 0.4 0.5 0.6 
Energy loss (eV) 



FIG. 1. Raman spectrum of YT1O3 at T = 13 K in (z, z) 
geometry, taken from Ref. |21] , A background is subtracted 
from the data. The sharp peak around 170 meV in the data 
is the two-phonon Raman signal, and is not considered in 
our theory. The thin-solid line is the superexchange theory 
curve. The anisotropy of the local model is reflected in its 
Raman spectra: (z, z) polarization (dashed line) gives a very 
different spectrum from (x,x) polarization (dotted line). In 
the superexchange model, the xx, yy- and zz-polarizations are 
equivalent. It should be noted that the experimental Raman 
spectra are also of cubic symmetry^. 



IV. RIXS 

The rapidly developing technique of Resonant Inelastic 
X-ray Scattering (RIXS) is an excellent probe of collec- 
tive excitations in transition metal compoundiPSJteSI The 
reason for this is that the X-rays carry enough momen- 
tum to map out the entire or at least a significant part of 
the Brillouin Zone, depending on the energy of the res- 
onant edge used. Due to the recent advances in energy 
resolution, it is now possible to access energies as low as 
~ 50 meV. This makes it in principle feasible to observe 
orbitons. 

In RIXS, the incoming X-rays excite a core electron 
into or above the valence band. Because it is a reso- 
nant technique, RIXS is clement-specific. Not only does 
this add more control to the experiments, it also helps 
with interpreting the RIXS spectra. In our case we con- 
sider the Ti 2p to 3d transition. This excitation can in 
principle affect the valence electrons in two ways: firstly 
through the core hole-excited electron pair's potential 
(from hereon referred to simply as core hole potential) 
and secondly, if the electron is excited into the valence 
band, by the Pauli exclusion principle. Because of these 
interactions with the valence electrons, the core hole can 
lose energy and momentum to the valence electrons. The 
core hole can affect the valence electrons on the core hole 
site itself, but it can also frustrate the bonds of this site 
with its neighbors. The intermediate state is shortlived, 
and when the photo-excited electron annihilates the core 



hole, the energy and momentum of the resulting X-ray 
photon are measured. From this measurement, it can 
be deduced what the energy and momentum are of the 
created excitations in the solid. 

In the experiment!^ we analyze, the L3 edge is used, 
where the 2p core electron is promoted from the spin- 
orbit split j = 3/2 state to a 3d state. The intermedi- 
ate states have a complicated multiplet structure, with 
large spin-orbit coupling in the core levels, strong intra- 
ionic Coulomb interactions altered by the core potential, 
etc, which makes the RIXS process hard to analyze mi- 
croscopically in an exact way. Fortunately, it is possi- 
ble to disentangle the problem of the intermediate states 
from the low-energy orbital transitions in the final states. 
Namely, since the intermediate states dynamics is much 
faster than that of orbital fluctuations, one can construct 
- based on pure symmetry grounds - a general RIXS 
operator describing orbital transitions between the ini- 
tial and final states. In this operator, the problem of 
the intermediate states can be cast in the form of phe- 
nomenological matrix elements that depend only on the 
energy of the incident photon and its polarization fac- 
tors. These martix elements can then be calculated in- 
dependently, e.g., by means of well developed quantum- 
chemistry methods on small clusters. This approach is 
general, but can be simplified in the (physically relevant) 
case where the energy dependence of martix elements is 
smooth: they can then be regarded as effective constants 
at energy scales corresponding to the low-frequency or- 
bital dynamics. 

RIXS spectra are described by the Kramers-Heisenberg 
formula, which can be written in terms of an effective 
scattering operator O q : 



A fi = (f\D 



Ea -H-iT 



:D\i) 



(/|O q j 



(29) 



where Ei is the incoming photon's energy, H is the 
Hamiltonian and D the dipolc transition operator. T is 
the lifetime broadening of the intermediate states. The 
cross section is obtained via 



do \ - 

dojdQ ^ 
f 



\A fi \ d(u-Wfi). 



(30) 



Here ujfi is the energy difference between the final and 
initial state of the solid. The cross section can also be 
written in terms of the Green's function for the effective 
scattering operator: 



do \ - 

dojdQ, ^ 

f 



J (/| 6 q \i) j' 6(u - ufi) = -hm {G(uj)} 

(31) 



with 



G(u) 



die*" (i\ 6 q (i)o q (0) |<) 



(32) 



The effective scattering operator can in general be ex- 
panded in the number of sites involved in the scattering 
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process: 



■Oi 



(33) 



where q is the transferee! momentum. The phase factor 
comes from the dipole operators. We neglect RIXS pro- 
cesses that create excitations on more than two sites in 
the final state, and further assume that the two-site pro- 
cesses are dominated by processes on nearest neighbors. 

One may distinguish two regimes for RIXS processes: 
in the first regime T is much larger than the relevant 
energy scales of the intermediate states, and these pro- 
cesses can be easily an alyzed with the Ultrashort Core 
hole Lifetime expansion^ 41 * 42 *. In the other regime V is 
small and its inverse is irrelevant as a cut-off time of the 
intermediate state dynamics. The lifetime broadening at 
the transition metal L edges is relatively small, and the 
effects of the core hole on the valence electrons is aver- 
aged over many precessions of the core hole due to the 
large spin-orbit coupling in the core levels of transition 
metal ions. Therefore the A\ g component of the Coulomb 
potential of the core hole dominates the scattering pro- 
cesses. In the following, we assume that the titanates 
belong to the regime of small V, and that the internal 
dynamics of intermediate states is the fastest process in 
the problem. 



Returning to the scattering operator, Eq. (33 1, we are 
left with two interesting cases. The single-site operator 
is dominated by the A\ g component, but this only gives 
contributions to the Bragg peaks. The subleading or- 
der therefore consists of single site processes Oi of other 
than A\ g symmetries, and of two-site processes Oy of 
A\ g symmetry. 

The single site coupling of RIXS to the orbitals can be 
dubbed a "shakeup" process. If we allow the core hole po- 
tential to have a symmetry other than A\ g , it can locally 
induce an orbital flip. If the orbital ground-state is dom- 
inated by superexchange many-body interactions, a local 
flipped orbital will strongly interact with the neighbor- 
ing sites and thus becomes a superposition of extended 
(multi-)orbitons. In the limit of strong crystal field split- 
tings, however, this excitation remains a localized, on-site 
transition between ti g levels. 

Two-site processes may involve modulation of the 
superexchange bonds, analogous to two-magnon RIXS, 
where the superexchan ge con stant J is effectively modi- 
fied at the core hole sit d 36 * 38 M ^. The core hole potential 
locally changes the Hubbard U, which in effect changes 
Jse = 4:t 2 /U on the Ti-Ti bonds coupled to the core 
hole site. Alternatively, the two-site processes can de- 
scribe the lattice-mediated interaction that is altered by 
the presence of a core hole. The equilibrium positions 
and vibration frequencies of the oxygens surrounding the 
core hole site may change, affecting the intersite inter- 
actions. As said above, the A± g component of the core 
hole potential is most relevant in the two-site coupling 
channel ? I . 



This section is divided into three subsections. Subsec- 
tion IV A deals with the single site shakeup mechanism 



and contains the evaluation in the superexchange model. 
The next subsection, |IVB| is devoted to the calculation 
of the same processes in the local model of the orbital 
excitations in YTi03 . The final subsection |IV C| covers 
two-site processes, evaluated within the superexchange 
model. A detailed comparison is made of the RIXS spec- 
tra arising from the different models. 



A. Single site processes — Superexchange model 

We start out with an analysis of the single site pro- 
cesses. RIXS processes that involve orbital excitations on 
a single site are dominated by direct transitions between 
the t2 g orbitals when the core hole potential is not of A\ g 
symmetry. In a superexchange dominated system, a lo- 
cal flipped orbital strongly interacts with the neighboring 
sites and becomes a superposition of extended orbitons. 

We start from the Kramers-Heisenberg equation 



0^ = D 



1 

Ej -H-iT 



D. 



(34) 



We insert the polarization-dependent dipole operator D 
which we take to be local: D = . Di with 



\m) (m\ r • e \d) (d\ 



d , m 

+ el qoufR 1 |^ ( d | f . £ > | m J ( m |) + h c ; (35) 



where e and e' are the in- and outgoing polarization vec- 
tors respectively, \d) denotes the state of atom i when it 
is not photo-excited and \m) denotes the system's inter- 
mediate eigenstates: 



H = E m H (m| 



(36) 



Now we consider only the single site part of the effective 
scattering operator in Eq. ( 33 1 : 



O t = \d')m'-e'\ 



1 



Ej — En 



iT 



e-r\d) (d\ 

Next we decompose the operator part into terms trans- 
forming according to the rows of the irreducible repre- 
sentations of the octahedral group (labeled by T, not to 
be confused with the core hole lifetime broadening): 



|d') (d\ = J2^'dr- 



(38) 



In the second quantized picture, we need only terms that 
are quadratic in the creation and annihilation operators. 
With the irreducible representations Ai g , T lu , E g and T 2g 
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all possible \d') (d\ can be constructed. Therefore f as- 
sumes only the following forms: 



M g 
T\ u 
E„ 



2;/ 



r = i 

fe{Q x ,Q z } 
f e{f x ,fy,f z }. 



(39) 
(40) 
(41) 
(42) 



The operators T and the corresponding 3x3 matrices T^^ 
are defined in AppendixlA"] Because A\ g only contributes 
to elastic scattering, we chop it from hereon. 

Further, we also decompose the dipole matrix elements 
into 



(43) 



{d'\P\m) (m\a\d)=J2 T ^ M *d 



with a, (3 <E {x, y, z } an d the M d , d listed in Appendix |b| 
Eqs. dBlJ through (JB9). Plugging Eqs. ([38]) and (|43j) into 
Eq. (37), we obtain 



d,d',m T' q,/3 1 m r 
which can be simplified using 

y"l M d'd T d'd = Srx' E Md'd^d'd- 



(44) 



(45) 



This identity can be proven by interpreting M r and T as 
matrices indexed by d and d'. Then it can be seen that 
M T <x T. We thus obtain 



M T d : d Y d , d = Tr (m t 't t ) ex Tr (rT 



(46) 



which is zero for r ^ f, proving the above identity. We 
find then 



6 i = y ^p r M v t i 

r 

with a polarization factor 

Pr = /J e^r/3Q£a 

a,B 



(47) 



(48) 



and the matrix elements Mr depending on the multiplct 
effects in the intermediate state 



d.d 



M d , d T d ' d 
? Ei-E m -iY 



(49) 



One can perform the sum over m, which yields 
M$ d 



E 



Ei — E m — iT 



and similar expressions for the other representations. As 
discussed above, we will assume that the intermediate 
state dynamics is much faster than that of t 2g orbitals 
we are interested in, and thus regard the matrix elements 
as phenomenological constants. Further, using that A^r 
does not depend on any coordinate and therefore must 
be invariant under the octahedral group, we obtain 



Mq = M Qx = M Qz 
Mt = Mt x = Mt v = 

Mi = Mr = Mi = Mr. 



(51) 
(52) 
(53) 



The Mr are hard to calculate explicitly since they involve 
inverting H, which contains the multiplct structure. In 
the following, we assume Mr = M for all V. This is a 
reasonable assumption: the core hole generates a multi- 
tude of many-body states that evolves very rapidly due 
to the large spin-orbit coupling and intra-ionic Coulomb 
interactions, and therefore its potential is averaged. Any 
particular symmetry is washed away; all become equal, 
except for the Ai g component, which is enhanced at the 
cost of the others. This is also the reason why the exper- 
iments at the t2 g and e g edges are similar^! the different 
edges create different multiplct structures initially, but 
these differences are averaged out by the intermediate 
state dynamics, as far as we are concerned with t 2g or- 
bital transitions at relatively low energies 0.2-0.3 eV. 

Note that Pr and Mr are independent of the site i. 
Only F depends on i, giving 



O q = ]TPrAlr]re"i- R -r\; 



(54) 



Most interference terms between different T's are zero. 
This comes about because of the specific ground state 
ordering. Transforming to the local axes (Eq. (15) in 
Ref. [5]), the ground state and Hq are invariant un- 
der translations, while the operators T a ^ and l a ^i (with 
a G {x, y, z}) acquire a phase upon translation to a differ- 
ent sublattice, which is equivalent to a momentum shift 
(by orbital ordering vectors) for the corresponding T q . 
Therefore, many interference terms are zero, which can 
be seen from Eqs. (31) and (32): two operators with 



different momenta cannot bring the ground state (zero 
momentum) back to itself. The only non-vanishing in- 
terference terms are (0| QJ. „(t)Q z „(Q) |0) which do not 

acquire momentum shifts and (0| q (i)/ a q (0) |0) where 
the momentum shifts cancel. 

To compare with experiment, we calculate the polar- 
ization factors Pr for the experimental setup of Ref. |23| 
where q is along the [001]-direction. Only the incom- 
ing polarization is fixed, the outgoing polarization is not 
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detected and should be averaged over. We have 

e = ( —= sin 9, —= sin 9, cos 9) 
V 2 v 2 

= ( -= sin 0, -= sin 0, cos 0) 

v2 v2 



(55) 
(56) 
(57) 



with 29 the scattering angle. Then, we find for the hor- 
izontal outgoing polarization e' H (i.e. the electric field 
vector is in the scattering plane): 



P 



Ptx,h — P' 

1 



P, 







Pa, h = - cos 29 



P 



1 



2V3 
1 . 



(1 + cos 2 9) 



Pt z ,h = -sm"V 
Pu.h = -Pl..h = 



V2 



sin a cos I 



(58) 
(59) 

(60) 

(61) 

(62) 



and for vertical outgoing polarization e' v (electric field 
vector perpendicular to the scattering plane): 



P Alg ,v = Pq z ,v = Pt z ,v = 
P Q*V = sin 6» 



Pt^v 
P*,v 



= -P 



1 



T V ,V 



2V2 

i 

2V2 



cos 9 
cos 9 



(63) 
(64) 

(65) 
(66) 
(67) 



For horizontal polarization, the polarization factors make 
all remaining interference terms zero. 

In Appendix O the one- and two-orbiton parts of 
the r q = J2i e are listed. They are obtained by 

performing the transformations on the orbital operators 
mentioned in Ref. [5]. Then, the %j> c orbital (with corre- 
sponding annihilation operator c) is condensed: 



n~ c = I c 1 + Sn s 



(68) 



where <5ng is the fluctuating part. In the completely or- 
dered state, (Sric) = (rig) = (rig) = and |co| 2 = 1, 
while in the completely disordered state (Sns) — {rig) — 
(n^) = 1/3 and |co| 2 = 0. Ref. [5] obtains a finite value 
for the quadrupole orbital order parameter: 



Q = n~ c -(ria + n h )/2 = {Qj + 5Q ~ 0.19 + SQ (69) 

with the fluctuating part averaging to zero. This fixes 
I Co 1 2 — 0.19. Taking the square root of Eq. (68 1, one 
arrives at 



to first order in the fluctuations Sn^. 

In the process of writing the T q in terms of orbiton op- 
erators, unphysical contributions to the intensity may ap- 
pear as a result of neglecting cubic and higher order terms 
in the orbiton operators. When restoring all terms, these 
unphysical contributions should cancel by symmetry. For 



q along the [001]-direction for instance, [H, Q z 



if 



we use the untransformed forms Eqs. ([IJ and (|A5j), and it 
is clear that there should only be an elastic contribution 
to the intensity. However, in terms of orbitons, this selec- 
tion rule is violated if we go only up to quadratic orbiton 
terms. To make sure these unphysical contributions are 
dropped, we first calculate the commutator in the un- 
transformed picture. If this yields zero, the commuting 
part of the scattering operator is dropped. Applying this 
procedure to the case where q is along the [001]-direction, 
we find that only Q z ,q among the operators (40 1 to (42 ) is 
zero while all the other channels give finite contributions. 

Since we did not include explicitly the orbiton-orbiton 
interactions, damping of the orbitons should still be taken 
care of, at least on a phenomenological level. As in the 
case of Raman scattering calculations, we introduce by 
hand an energy broadening 7 of the orbiton states (half- 
width at half maximum, HWHM) of 7 = 0.4 J orb . This 
broadening can also be used to take orbiton damping by 
phonons, magnons etc. into account. In addition to this, 
there is an experimental broadening added of 27.5 meV 
(HWHM)P. 

The resulting spectra are shown in Fig. [2j The 
intensity is strongly momentum-dependent (especially 
for q along the z-direction) , which is also seen in the 
experiments^!. This dependence is mainly due to the co- 
herent response of the exchange-coupled orbit als which 
enhances at large momenta, reflecting staggered orbital 
order in the ground state - Eq. ([3|. In Fig. [3] the the- 
oretical cross section (with q along the z-direction and 
horizontal incoming polarization, i.e. the electric field is 
in the scattering plane) is compared to the experimental 
data. The main features of the data^ are reproduced: 
the spectral weight increases with increasing q z and there 
is virtually no dispersion of the maximum of the theoret- 
ical curve (because it is determined by the two-orbiton 
continuum, containing an integration over the Brillouin 
zone) . 

Especially in the second and third plots, the one- 
orbiton shoulder seems a bit too large. However, we note 
that there are several factors that can alter the line shape. 
First we note again that the weight of this shoulder is con- 
trolled by the orbital order parameter: if the orbital order 
melts, I Co 1 2 decreases and the one-orbiton peak becomes 
less intense. The value we used (|c | 2 = 0.19) is obtained 
at zero temperatur e, a ssuming that YTiC>3 is a fully sat- 
urated ferromagnetP^l Under realistic conditions, |co| 2 is 
expected to be smaller than 0.19. Indeed, the saturated 
ordered moment in YTiC>3 is actually 0.84/i^, which is re- 
duced further to approximately 0.80/is at T — 15 Kl 27 l 43 [ 
Correspondingly, the orbital order is decreased by joint 
spin-orbital quantum fluctuations, suppressing the one- 



5 f = V M 2 + Sn E « |c | + — Sn- C (70) 
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orbiton peak. 

Secondly, we assumed all A^r are equal. Different val- 
ues would correspond to different line shapes. We note 
that the Ti g representation has a much reduced one- 
orbiton contribution compared to the other channels. 

Thirdly, we introduced the finite orbiton lifetime 
broadening as a phenomenological damping only. All ver- 
tex corrections to the two-orbiton diagram are neglected. 
In analogy to two-magnons, these terms can give correc- 
tions to the spectrum. 

The best chance to see a one-orbiton contribution to 
the spectrum is with momentum transfer directed maxi- 
mally in the [110]-direction. Fig. [4]shows the prediction 
for the shakeup mechanism with |co| 2 = 0.19: the one- 
orbiton peak is about as strong as the two-orbiton peak. 



B. Single site processes — Local model 

Although the response functions of the local model 
of YTi03 are entirely different from the superex- 
change model, the phenomenological scattering operator 
Eq. (|33| is still valid. Focusing on single site processes, 
Eq. ( |54[ ) can be evaluated using the wave functions found 
by Pavarini et al.: Eqs. (|5| through ([7]). Since the eigen- 
states of the local model have a very simple form, we can 
straightforwardly use Eqs. ( |A1| through ( |A8| to evaluate 
the RIXS spectrum. The Pr and Afr remain the same 
as in the collective orbiton case. The spectrum now con- 
sists of two sharp peaks at u>± and u>2- These peaks can 
be broadened by coupling to the lattice as well as due to 
the superexchange coupling. 

Because there are four sublattices which all support 
their own, local eigenstates, the RIXS intensity can be 
deco mpo sed into four signals. From the expressions 
(12 1 through (15), it is easily derived how the Tj 



Eqs. 
transform. 

So far, the analysis is similar to Sec. IV A However, in 
the local model, the eigenstates are local and this changes 
the analysis of Sec. |IV A| at two important points. The 
first one is that the momentum shifts of T q do not destroy 
interference terms: any final state |/) can be reached with 
any shift of q. All interference terms can in principle be 
present. The second point to be noted is that, because 
the eigenstates are local, the only momentum dependence 
of the cross section comes in through the experimental 
geometry, which is reflected in the polarization factors 
Pr- 

When we compare the theoretical RIXS spectrum of 
this model to experiment, we again have to take into ac- 
count the average over the two outgoing polarizations. 
Assuming again Air are the same for all T, and intro- 
ducing the same broadening as before (HWHM 7 « 30 
meV phenomenological intrinsic broadening plus 27.5 
meV HWHM experimental broadening), we obtain the 
spectra shown by the solid lines in Fig. [5] 

It is evident that the local model yields a RIXS spec- 
trum that does not agree well with experiment. Firstly, 
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FIG. 2. (Color online) RIXS spectra for a model of 
superexchange-driven orbital order with RIXS coupling to 
orbitons only via the single site mechanism. The first two 
spectra are for q directed along the [001]-direction, the last 
two for q along the [110]-direction. The first and third spec- 
trum are for horizontal incoming polarization (electric field in 
the scattering plane), the second and fourth are for vertical 
incoming polarization (electric field perpendicular to the scat- 
tering plane) . Note that the q = points are different in each 
spectrum because of the different experimental geometries, 
leading to different Pr- We only plotted the experimentally 
accessible part of the Brillouin Zone. 
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q z = 0.24 Jt q = (0.36 Jt,0.36 Jt,0) 

H-pol in 

V-pol in 




-7 -6 -5 -4 -3 -2-10 1 -7 -6 -5 -4 -3 -2 -1 

Energy loss (G0/J orb ) Energy loss (oo/J orb ) 



q z = 0.40 Jt 
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q z = 0.52 Jt 
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FIG. 3. RIXS spectra for a superexchange-driven orbital 
order with RIXS coupling to orbitons only via the single 
site mechanism (solid line), compared to the experimental 
datd^SI. The vector q is directed along the [001]-direction, 
with q z as indicated in the figures. We took J or b = 80 meV, 
and introduced a phenomenological HWHM broadening of 
7 = 0.4 J or b ~ 30 meV for the orbitons, as well as the HWHM 
experimental resolution of 27.5 meV. The elastic peak is fitted 
with a Gaussian (dash-dotted line). 



FIG. 4. Spectra, obtained with the single site mechanism, 
for the largest experimentally accessible momentum transfer 
directed along the [110]-direction. The solid line indicates 
the case where the incoming polarization is horizontal, the 
dashed line is for vertical incoming polarization. The elastic 
peak has been removed. In the horizontally polarized case, 
the single orbiton peak is quite strong and should be visible 
in experiments if the system is superexchange-driven and the 
RIXS signal is dominated by the single site mechanism. 



there is no two-peak structure visible in the data. The 
presence of a two-peak structure in the theoretical curves 
does not depend on the assumption that all the .Mr are 
equal. We may finetune the model to produce a better 
fit by changing the energy levels found in Ref. [11] so that 
both crystal field transitions have an energy of 240 meV, 
and introducing a very large intrinsic broadening of 100 
meV (see the dashed lines in Fig. [5]). But even in the ar- 
tificially optimized case of degenerate levels to produce a 
single peak, the intensity trend remains in contradiction 
with experiment. Further, it is impossible to tune the 
energy levels to optimize simultaneously the RIXS and 
Raman data. Both experiments show a peak at the same 
energy, while the local model theory predicts the Raman 
spectra (with its double crystal field excitations) to peak 
at approximately double the RIXS peak energy. 

Even though we could improve the line shape by in- 
creasing 7, the intensity gain with increasing q z cannot 
be reproduced in any way. In fact, the trend is the oppo- 
site: as q z increases, the spectral weight of the theoret- 
ical spectrum decreases (see Fig. [5]) . We recall that the 
q-dependence in this case is merely due to polarization 
factors Eqs. (58-671, since in a local picture, each Ti ion 

This is 



contributes independently to the cross section, 
in sharp contrast with the superexchange picture, where 
the intensity has an intrinsic q-dependence because of 
the collective response of all the Ti ions. 

Finally, it is hard to reconcile the temperature depen- 
dence of the experimental data with the local model. The 
peak is seen to broaden and lose a large part of its spec- 
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FIG. 5. RIXS spectra for the local model (solid and dashed 
lines) compared to experimental datcP^. The vector q is di- 
rected along the [001]-direction, where q z is indicated in the 
figures. The dashed curve shows the artificially optimized 
model with degenerate crystal field levels. We introduced a 
phenomenological intrinsic HWHM broadening (7 « 30 meV 
[solid line] and 7 = 100 meV [dashed line]) for the final states 
and added experimental broadening. 



tral weight with increasing temperature^. Ascribing this 
broadening to phonons has two difficulties. Jahn- Teller 
active phonons have energies around 30 — 60 meV and 
are therefore not very sensitive to temperature up to 
T s» 350 — 700 K. Further, such a broadening would im- 
ply a strong orbital-lattice coupling. This then raises the 
question why no structural phase transition is seen in the 
titanates. We note that this is very different from, e.g., 
manganites where orbital-lattice coupling dominates. 



C. Two-site processes 

The second term in the expansion of the effective scat- 
tering operator, Eq. ( |33[ > , involves two-site processes. 
Due to the strong multiplet effects, the core hole poten- 
tial is averaged out and becomes mainly of A\ g symmetry. 
While such a potential cannot directly flip the orbitals 
at the core hole site, it does affect multi-site processes. 
In the case of the superexchange model, the core hole 
potential effectively changes the superexchange constant 
Jse locally as disc ussed earlier in the context of RIXS 
on magnon d 36 * 38 * 40 ! This process is illustrated in Fig. [6] 
(In principle, it is also possible that the core hole poten- 
tial modifies the orbital interactions via the lattice vi- 
brations.) In this section, we consider the superexchange 
modulation mechanism to illustrate two-site process in 
RIXS on orbital fluctuations. 



time 



tv 

-o- -< 
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3d to. 
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2p -O- 

Ti 



Ti 



FIG. 6. (Color online) In the superexchange model, two-site 
RIXS processes locally modify the superexchange interaction, 
coupling the RIXS core hole to the ti g orbitals. Shown is 
an orbital superexchange process between two neighboring Ti 
ions in the presence of a core hole. On the left, one of the ions 
is excited by an incoming x-ray photon. After that, the t^ g 
electrons undergo a superexchange process. On the right, the 
virtual state of the superexchange process is depicted. The 
presence of the core hole frustrates the superexchange process. 
Instead of the usual Hubbard U, the energy of the virtual state 
is lowered by the presence of the positively charged core hole. 
This modifies the superexchange constant Jse — 4t 2 fU at 
the core hole site. 



The superexchange modification can be derived explic- 
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itly by starting from a Hubbard model 

i 

- Uc^PiPj (fW + n-b,i + n c>i - 1) (71) 

where the last term includes the Coulomb energy U c of 
the core hole attracting the t 2g electrons, pi is the anni- 
hilation operator for 2p core electrons at site i. We have 
taken the core hole potential to be of A\ g symmetry. Do- 
ing perturbation theory to second order in t/U( c ) (U and 
U c are about the same order of magnitude), we obtain 
the superexchange Hamiltonian 

H = H +J2PiPl UA^ +S - + const.) (72) 



i,6 



with <5 pointing to nearest neighbors, nf^ — n a ^ + nt,^ 
(the other can be obtained by permuting the indices 
a, 6, c) and 



Ji = 
J 2 = 



t z 



u -u c u 



u + u c u-u c 



2^ 

IT 



(73) 
(74) 



so that Ji/ J 2 = (1 + J7/f/ c ) /2- Eq. ( 72 1 shows we get the 



unperturbed Hamiltonian plus a contribution which is 
active only if there is a core hole (in which case Pip\ — > 1). 
The J\ term involves single site processes only, and is 
therefore included in the general description in Sec. |IV A] 
In the following, the J\ term will be dropped. 

For simplicity, polarization effects arc neglected and 
we assume U c to be independent of the specific dipolc 
transition. We take 



E U'^^Pi 4 + e^^p\d t ) + h.c. (75) 



D 



with pi the 2p electron annihilation operator and di the 
3d e g electron annihilation operator. The position of the 
i th site is Hi. The transfered momentum is q = q ou t— qin- 
The neglected polarization dependence could give rise to 
a q-dependent factor in the cross section, but will not 
affect the line shape for a specific q. 

The relevant energy scale for the excitation of orbitons 
in the intermediate states via superexchange bond mod- 
ulation is J 2 , as established above, as long as the core 
hole potential is of A\ g symmetry. This is the case when 
core level spin-orbit coupling and Hund's rule coupling 
are large compared to J 2 : the core hole evolves rapidly 
with time and its potential's symmetry averages out to 
A\ g before any orbitons can be excited. Because the sym- 
metry is effectively cubic, bonds in all directions are af- 
fected in the same way. The effective scattering operator 



must therefore be a function of J^s -^-ii+s which is of A\ g 

symmetry. Non-linear operators like A^ +S A^ +S , are ex- 
cluded, they are expected to yield smaller contributions 
because more and more distant sites are involved. In 
the expansion Eq. (33 1, these come in at different orders. 



The only remaining candidate for the two-site effective 
scattering operator is therefore 



O a = M, 



i,6 



D iq-R iy |(7) 



i+S 



(76) 



where M. 2 is an unknown phcnomcnological matrix ele- 
ment, in the same way as in Sec. |IV A| By construction, 
the two-site process matrix element Ai 2 should be pro- 
portional to J 2 with a constant determined by the inter- 
mediate state dynamical susceptibilities. At this stage, 
without microscopical calculations of the single-site .Mr 
(49 1 and two-site M. 2 matrix elements, we cannot judge 



which coupling process dominates the observed RIXS on 
orbital excitations. Instead, we calculate two-site pro- 
cess independently and compare it with both experimen- 
tal data and the results obtained above for single-site 
coupling mechanism. 

As it turns out, the two-site effective scattering opera- 
tor ( 76 ) contains only two-orbiton creation terms; it does 



not create single orbitons because the orbitons are con- 
structed in the first place to diagonalize the Hamiltonian: 
all linear contributions to Aff in Eq. (j^) are canceled 
(similar to the Raman scattering calculations above). 

Using again the transformations on the orbital oper- 
ators mentioned in Ref. |5J, condensing the ip c orbital 
and transforming to orbiton operators, we obtain for the 
two-orbiton creation part 



^ 2) =-^ 2 E[Ai(k,q)4k«t- 



k-q 



/ 22 (k, q)4,k4,-k- q + ii2( k > q)"i,k a 2,-k 



t „,t 



(77) 



where the (k, q) are lengthy functions listed in Ap- 
pendix ID] The cross section then is 



' ! ~ a> " } oc E I </l O q 2) |0) J 2 S(oj - w 1/2 , k - w 1/2)k+q ) 



= i^[|/ 11 (k,q)+/ 11 (-k-q,q)| 2 

k 

5(<J - Wl,k - ^l.k+q) 

+ l/22(k,q) + / 2 2(-k-q,q)| 2 x 

S(0J - W 2 ,k - ^2,k+q) 

+2 |/i2(k, q)| 2 8{u) - Wi, k - ^2,k+ q ) 



X 



(78) 



The resulting cross section for transfered momenta along 
the [001] direction is shown in Fig. [7] As in the above sec- 
tions, we introduced here by hand an energy broadening 
7 of the orbiton states of 7 = 0.4 J or t,. 
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(0,0,Jt) 



FIG. 7. (Color online) RIXS spectrum for two-site processes 
within the superexchange model. The color scale denotes 
the intensity. The figure shows q running from (0, 0, 0) to 
(0, 0, 7r). One-orbiton creation is not allowed for the superex- 
change modulation mechanism. The intrinsic energy broad- 
ening 7 of the orbiton states is 7 = 0.4 J or b » 30 meV, and 
the added experimental resolution of 27.5 is approxi- 

mately 0.34 Jorb. 



A few things should be noted. Firstly, the spectrum 
disappears at q = 0. This is clear from Eq. (76 1: the 



scattering operator becomes proportional to the Hamil- 
tonian Eq. (JlJ, giving elastic scattering only. 

Secondly, the spectrum shown in Fig. [7] is calculated 
without taking polarization dependence into account. 
That could change the relative spectral weight for dif- 
ferent q's, but does not affect the line shapes. 

In Fig. [8] we compare the calculated superexchange 
spectra for the specific q values of the experiments re- 
ported in Ref. [23]. The only free parameter (J or h) gives 
a best fit for J or h = 75 meV. As is evident, the increase in 
spectral weight is qualitatively accounted for by the the- 
ory, although the theoretical curves show a much stronger 
increase with increasing q z . We note that one factor that 
could diminish this discrepancy is, as stated above, the 
polarization factor we omitted: it could change the rela- 
tive weight (but not the line shape). 

Further, when we compare the theoretical line shapes 
with the experimental ones, the high energy tail of the 
experimental data is a bit more intense than in our calcu- 
lations. This could perhaps be accounted for if one would 
consider multi-orbiton scattering. Likewise, the shoulder 
around u> — 1.5J or h ~ 110 meV could be due to two- 
phonon processes. Both these discrepancies depend on 
the choice of 7: a larger orbiton damping would transfer 
spectral weight from the center of the theoretical peak to 
its tails. 

Summarizing, two-site processes can capture some of 
the features seen in the RIXS data (the intensity trend 
with increasing momentum transfer, and a single peak 



without dispersion) , but the overall fit is less satisfactory 
compared to the results of the single-site process shown 
in Fig. |3J 



V. CONCLUSIONS 

We have considered two different models, widely dis- 
cussed in literature to describe orbital physics in titan- 
ites, in the context of Raman and x-ray scattering exper- 
iments. These models correspond to two limiting cases 
where the orbital ground state is dominated either by col- 
lective superexchange interactions among orbitals or by 
their coupling to lattice distortions. The models predic- 
tions, obtained within the same level of approximations, 
are compared to the experimental data on Raman (Fig.l) 
and on x-ray (Figs. 3 and 5) scattering in titanites. What 
is evident from this comparison and our detailed analy- 
sis is that the local crystal field model of YTiC>3 fails to 
give a coherent explanation of both Raman and RIXS 
data taken together. There is no way one can get rid of 
the two-peak structure predicted for RIXS by this model 
without artificially finetuning its parameters. Further, 
once tuned to the RIXS spectra, the Raman spectra 
will be impossible to fit with the local model anyway, 
since it yields double dd-excitations, different from the 
single crystal field excitations in RIXS. Experimentally, 
however, both techniques show a peak at the same en- 
ergy. Also, a huge anisotropy between out-of-plane and 
in-plane polarizations is predicted by the local model, 
which is not observed in Raman data. Further, the tem- 
perature dependence of the experimental data is hard to 
explain from a local model: the intensity of crystal field 
transitions is expected to remain unchanged. Finally, the 
q-dependence of the RlXS-intensity is not reproduced by 
the local model; in fact, the trend is opposite. We be- 
lieve especially the last four points are robust evidence 
that the 250 meV peak seen in Raman and RIXS is not 
due to local dd-excitations. 

On the other side, the picture of collective excitations 
offers much better and broad agreement with the experi- 
mental data. The general features of both the Raman 
and RIXS data are reproduced by the superexchange 
model. For RIXS we presented a phenomenological scat- 
tering operator for single and two-site processes, evalu- 
ated within the superexchange model. Although both 
single and two-site processes the general trends of the 
RIXS data right, the two-site processes clearly have a too 
strong q-dependence of the intensity. The RIXS spectra 
obtained with the single site operator fit the data very 
well, suggesting that this process of generating orbitons 
might be the predominant one in the transition metal 
oxides. The only slight deviation from the experiments 
is the one-orbiton peak, which our theory overestimates. 
However, we note that the theoretical one-orbiton peak is 
decreased if we consider realistic conditions such as finite 
temperature and residual spin-orbital quantum fluctua- 
tions, which obstruct the orbital order and reduce |cq| 2 , 
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FIG. 8. Theoretical RIXS spectra for two-site processes, cal- 
culated within the supcrexchange model (solid line), com- 
pared to experiment!^, q is directed along the [001]-direction, 
with q z as indicated. We obtain a best fit for J or b = 75 meV. 
The solid lines are cuts from the plot of Fig. [7] where we added 
a Gaussian fit to the elastic peak (dash-dotted line). A phe- 
nomenological intrinsic orbiton broadening of 0.4 riJsE ~ 30 
meV is added, as well as the experimental resolution of 27.5 
meV (both HWHM^H. 



which in its turn controls the one-orbiton spectral weight. 

Ref. [23 reported that the RlXS-intensities in LaTiO,3 
and YTiC>3 show different q-dependences: While the 
intensity in YTi03 increases with q, it decreases in 
LaTiOa . On a qualitative level, this contrasting behavior 
can be understood from the superexchange picture as a 
manifestation of the (dynamical) Goodenough-Kanamori 
rules, according to which the spin and orbital correlations 
are complementary to each other. This implies that the 
spin and orbital susceptibilities are expected to behave in 
an opposite fashion. Since magnetic orderings in YTi03 
and LaTiOs are different (ferro- and antifcrromagnetic, 
respectively) , collective response of orbitals in these cou- 
pounds are expected to be enhanced also at different por- 
tions of the Brillouin zone: at large q in YTi03 and, in 
contrast, at small q in LaTiOs, which are complemen- 
tary to the respective locations of their magnetic Bragg 
peaks. The superexchange picture suggests also that the 
q-dependence of the orbiton RIXS- intensity should have 
cubic symmetry in both LaTiOjj and YTiC>3, as follows 
from their isotropic spin-wav e! 6 * 35 ^ and Raman spectra^. 
Future RIXS experiments in titanites would be useful to 
verify these expectations. 

A previous estimate^ from neutron spin wave data^ 
puts the orbital exchange constant J or b at 60 meV. In 
close agreement with this estimate, the theoretical Ra- 
man spectrum fits best to experiment when J or b = 
65 meV (vertex corrections may change this number, 
though). Matching to a lesser degree to the estimate, 
we find for the RIXS spectra J orb = 75 and 80 meV for 
the two-site and single site processes, respectively. 

To establish the nature of the 250 meV peak, it is 
of great importance to search for the one-orbiton peak. 
In Raman scattering one-orbiton creation seems to be 
strongly suppressed, but in RIXS it would be possible 
to see a one-orbiton peak when q is directed maximally 
along the [110]-direction. There the one-orbiton peak 
(around u ~ 1.8 J or b ~ 140 meV) is about as strong 
as the two-orbiton continuum, assuming single site pro- 
cesses are the dominant RIXS channel, and provided |co| 2 
is not too small. 

To summarize, we may conclude that the existing Ra- 
man and RIXS data in titanites are better described by 
the superexchange model. This implies that while some 
polarization of orbitals by static lattice distortions must 
be a part of a realistic, "ultimate" model for titanites, 
the orbital fluctuations which are intrinsic to the t^g or- 
bital superexchange process^ are not yet suppressed and 
strong enough to stabilize nearly isotropic charge distri- 
butions around the Ti-ions. 

On the technical side, we believe that our semi- 
phenomcnological approach to the RIXS problem which 
disentangles the high-energy intermediate state dynam- 
ics from low-energy collective excitations of orbitals and 
spins may serve as a simple and efficient tool in the theo- 
retical description of Resonant Inelastic X-ray Scattering 
in oxides in general. 
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Appendix A: RIXS - Single site processes 

The angul ar mom entum I and quadrupole operators 
Q,f in Eqs. (40-421 are defined as follows: 



l x = i(c<b- & f c) 
l y = i(a/c — c'a) 
l z = i{b^ a — a'b) 

rib - n a 



1 



()_=~(P X + P y -2P Z ) 



(2n c 



T 



VyV z 



l Z ly 



V3 
{tfc + Jb) 

T y = l x l z + l z l x = -(c f a + a f c) 

-(a^b+tfa) 



rib) 



T, 



Ixly 



lyl X 



(Al) 
(A2) 
(A3) 
(A4) 

(A5) 

(A6) 
(A7) 
(A8) 



which are normalized by Tr (f 2 ) —2. The correspond- 



ing matrices IV d in Eq. (38 1 are 

T T * = - 
T T *=- 




i 0> 

-i | (A9) 
0, 



with the indices d, d' = (yz, zx, xy) (or for polarization 
dependence: a,0 = (x, y,z)). 

Appendix B: Multiplet factors 



For the multiplet effect factors in Eq. ( 43 1 , we have 



Kd 9 = Jl ((d'\ x \m) (m| x \d) + (d'\ y \m) (m\ y \ 



■ (d | z \m) (m| z \d)) 



(Bl) 



M%* d = ((d'\ y \m) (m| y \d) - (d'\ x \m) (m\ x \d)) (B2) 



— (2 (d'\ z \m) (m| z \d) - (d'\ x \m) (m\ x \d) 
v3 

-(d'\y\m)(m\y\d)) (B3) 
Mjr d = - ((d'\ y \m) (m\ z \d) + (d'\ z \m) (m| y \d)) 



(B4) 



M^ d = - ((d'\ z \m) (m\ x \d) + (d'\ x \m) (m\ z \d}) 



(B5) 



M ld = - i(d'\ * H (m| y \d) + (d'\ y |m) (m| x \d)) 



(B6) 

M\l d = -i({d'\ y \m) (m\ z \d) - (d'\ z \m) (m\ y \d}) 

(B7) 

M l d v , d = -i({d'\ z \m) (m\ x \d) - (d'\ x \m) (m\ z \d)) 

(B8) 

M l d 1 d = -i((d'\x\m) (m\y\d) - (d'\y\m) (m\x\d)) 

(B9) 

Note that the position operators act on the core electrons, 
not the tig ones. Both the core and ti g electrons are 
implied in the states \d) , \d'). 



Appendix C: The operators f in terms of orbitons 

In terms of the orbiton operators, we obtain the one- 
orbiton creation part of T q = Y^i e lq ' ri Fi to be 



1G 



= [{(! - ^Ki+q + (1 + V3)« qi+q } (sh 0i, qi+q + ch 0! 

+ {(-1 - \/3)w qi+q + (1 - \/3)w qi+q } (sh e 2 

,qi+q + C h #2, qi + q )a 2 

^ q = ^ a/ f [{(1 + \/3N 2+q + (1 - \/3) Uq2+q } (sh 1>q2+q + ch 0i, q2+q )c4,- q2 - q 



qi+q;"l,-qi-q 
2,-qi-q 
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2,-q 3 -q 
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{(1 - V3)uq 1+q + (1 + \/3> qi+q } (ch 2 

,qi+q _ S h 02, qi +q)ai2, 



'2,-qi-q 



j-(i) = ir^ii^/Tv 
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1 " ^ r {(l - \/3>q 2+ q + (-1 - \/3)« q2+q } (ch liq2+q - sh 9 lM2+cl )al 

|(1 + \/3)u q2+q + (1 - \/3)Wq 2+q | (Ch 2 , q2+q - sh 2 ,q 2+q )c4 _ 



-q2-q 



co 



z,q 
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qa-q 



Qi 1 
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with qi = (7r, 0, 7r), q 2 = (n, ir, 0), q3 = (0, ir, ir). The expressions for the two-orbiton creation part of f q = J2i e' q r *f , 



are 



'l 2 q = ^ X - uv')ch 0i sh ffi aj k aj _ k _ qi _ q + (vu' - uv')ch 2 sh 2 4,k4 ,-k- 



+ (mt' + w/)(ch 0! sh 2 - sh 0i ch fl^i.k^.-k 



-qi-q 



(C9) 



with u,v,9 1 ,6 2 = u k ,Wk,^i,k,^2,k and u',v',6[,9 2 = u k + qi +q, Uk+ qi +q, 0i,k+ qi +q, 02,k+ qi +q- Further, ^ 2 q and fi 2 q 
have the same form as fi 2 q but with qi replaced by q 2 and q 3 respectively. Next, 
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where in the expression for Tj£ we replaced qi by q 2 : U , v, 0[, 9' 2 = u k+q2+q , u k+q2+q , 0i ik+q2+q , 2 ,k+q 2 + q - 



f *« = E + vv ') - s {uv ' + W) ) ch 9i sh ^ Q tk«l-k- q3 - q 

+ j-(W + TO') + ^(W + Vu')\ Ch 6» 2 Sh 9' 2 a2,k a 2,-k-q 3 -q 

+ |-(W - VU') + ~(W - TO') | (ch 0! Sh 6> 2 + Sh 9 X Ch ^2) a l,k a 2,-k-q 3 -q 

where we replaced qi by q 3 : u' , v', 0[, 9' 2 = u k+q3+q , v k+q3+q , 6> ljk+q3+q , #2,k+q 3 + q - Finally 



(2) 



Qt 
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where in both equations we replaced qi by 0: u',v', 9[, 9 2 = Mk+g , ^k+g , fli,k+ q , 02,k+g • 
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Appendix D: RIXS - two-site processes with 
superexchange model 



Functions fn,f 22 and /i 2 in Eqs. (77-781 are: 



/n(k, q) = [-73, q (W + u'v) - 72, q (uu' - to') - (1 + 7i, q )(«u' + to')] (ch 9 X sh 9\ + sh 9 X ch 9\) 

+ 2 ["/[(uu' + to') + 7 2 (m/ - to') + j' 3 (uv' + u'v)} (sh X sh 9\ + ch 6>i ch 9[) (Dl) 

/2 2 (k, q) = [j3,q(uv' + u'v) + 72, q (W - w') - (1 + 7i, q )0"' + «">')] (ch #2 sh 9' 2 + sh 6> 2 ch 9' 2 ) 

+ 2 [7i(W + to') - 7 2 (W - to') - 73 (uv' + u'v)} (sh 9 2 sh 6> 2 + ch 9 2 ch <9 2 ) (D2) 

/i 2 (k, q) = 2 [73, q (W - uv') + 72,qH' + u'v) - (1 + 7i, q )(™' - u'v)] (ch 6»i sh 9' 2 + sh 6>i ch 9' 2 ) 

+ 4 - u'v) + j' 2 (uv' + u'v) - j' 3 (uu' - to')] (sh 9 Y sh 9' 2 + ch Q x ch <9 2 ) (D3) 

where we shortened notation by writing 9 1/2 = 6» 1/2 , k , = #1/2, k+ q , "(') = «k(+q), «(') = u k(+ q ) and 7- = 7 llk+q . 
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